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RESUMEN 

Calculamos la Zona dc Habitabilidad Galactica (ZHG) dc M31 basandonos, 
principalmente, en la probabilidad de formation de planetas terrestres depen- 
diente de la mctalicidad (Z) del medio interestclar. La ZHG fuc determinada 
a partir de un modelo de evolution qm'mica construido para reproducir un gra- 
dicnte dc mctalicidad prcciso en el disco galactico: [0/H](r) = —0.015 ±0.003 
dex kpc -1 x r(kpc) +0.44 ± 0.04 dcx. Estc gradicnte es el mas probable en 
presencia de dispersion intrinseca en los datos observacionales. El modelo de 
evolution quimica prcdicc una formation cstelar mas activa y una formation 
galactica menos efitiente, comparada con las de la Via Lactea. Suponicndo 
que los planetas tipo Ticrra sc forman bajo una ley de probabilidad que siguc 
la distribution dc Z mostrada por las estrellas con planetas, entonces la ZHG 
mas probable con vida simple se localiza entre 6 y 17 kpc en planetas con 
edades entre 1 y 4.5 Giganos, mientras que la ZHG mas probable con vida 
complcja sobreviviente a explosiones de supernovas se localiza entre 3 y 13 kpc 
en planetas con edades entre 4.5 y 6.5 Giganos. Dcbido a que la historia de 
formation estelar de M31 ha sido mas cfiticntc y que su ZHG en mas amplia 
y vieja que en nuestra galaxia, el mimero de planetas que puede albergar vida 
en M31 podrfa ser mayor que en nuestra Galaxia. 

ABSTRACT 

We have computed the Galactic Habitable Zones (GHZs) of the Andromeda 
galaxy (M31) based mainly, but not exclusively, on the probability of terres- 
trial planet formation, which depends on the metallicity (Z) of the interstellar 
medium. The GHZ was therefore obtained from a chemical evolution model 
built to reproduce a precise metallicity gradient in the galactic disk, [0/H](r) 
= -0.015 ± 0.003 dex kpc" 1 x r(kpc) +0.44 ± 0.04 dex. This gradient is the 
most probable when intrinsic scatter is present in the observational data. The 
chemical evolution model predicted a higher star formation history in both 
the halo and disk components of M31 and a less efficient inside-out galactic 
formation, compared to those of the Milky Way. If we assumed that Earth- like 
planets form with a probability law that follows the Z distribution shown by 
stars with detected planets, the most probable GHZ with basic life is located 
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between 6 and 17 kpc on planets with ages between 4.5 and 1 Gy, and the 
most probable GHZ with complex life that survived supernova explosions is 
located between 3 and 13 kpc on planets with ages between 6.5 and 4.5 Gy. 
Since the star formation history of M31 has been more efficient and its GHZ 
is wider and older compared to those of the Milky Way, the number of planets 
harboring life in M31 may be higher than that of our galaxy. 

Key Words: Galactic habitability — Chemical evolution — Abundance gra- 
dients — Andromeda galaxy — M31 

1. INTRODUCTION 

The Galactic Habitable Zone (GHZ) is defined as the region with sufficient 
abundance of chemical elements to form planetary systems in which Earth-like 
planets could be found and might be capable of sustaining life (Gonzalez et 
al. 2001, Lineweaver 2001). An Earth-like planet is a rocky planet character- 
ized in general terms by the presence of water and an atmosphere (Segura & 
Kaltenegger 2009). 

GHZ research has focused mainly on our galaxy, the Milky Way (MW) 
(Gonzalez et al. 2001, Lineweaver ct al. 2004, Prantzos 2008, Gowanlock et 
al. 2011). Gonzalez et al. were the first to propose the concept of a Galactic 
Habitable Zone, which is a ring located in the thin disk that migrates outwards 
with time. Lineweaver et al. (2004) later proposed that the Milky Way's GHZ 
is a ring located in the Galactic disk within a radius interval of 7 to 9 kpc 
from the center of the MW, and that the area of the ring increases with the 
age of the Galaxy. On the other hand, Prantzos concluded that the current 
GHZ is a wider ring that covers practically the entire MW disk. These studies 
confirmed that the Solar System is located within the GHZ since the Sun is 
found at 8 kpc from the center of the Galaxy. Nevertheless, a star with an 
Earth-like planet capable of sustaining life is more likely to be found in inner 
rings of the Galactic disk, between 2 and 4 kpc, owing to the high stellar 
surface density that presents the inner disk of our galaxy (Prantzos 2008, 
Gowanlock et al. 2011). 

The GHZ has recently been computed for two elliptical galaxies (Suthar 
& McKay 2012). Imposing only metallicity restrictions for planet formation, 
these authors found that both elliptical galaxies could sustain broad GHZs. 

Here, we extend MW studies to the disk of the most massive galaxy in the 
Local Group: the Andromeda galaxy (M31). M31 is a type SAb spiral galaxy, 
whose visible mass is 1. times larger than that of the MW, and M31 is at a 
distance of 783 ± 30 kpc (Holland 1998). 

The GHZ depends mainly on the abundance of chemical elements heavier 
than He (metallicity, Z), since Z leads to planetary formation. Moreover, the 
GHZ also depends on the occurrence of supernovae (SNe), since SNe determine 
regions where planets can survive possible sterilization. 

In this paper, we excluded the bulge of Andromeda from the GHZ, despite 
its having a high enough abundance of chemical elements to form planets 
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(Sarajedini & Jablonka 2005). The bulge does not provide a stable habitat 
for life due to its high supernovae rate, which could sterilize planets (Meneses- 
Goytia 2009), and the proximity of stars that can destabilize the orbits of the 
planets (Jimenez-Torres et al. 2011). 

Throughout this paper the terms "evolved life" and "complex life" are 
synonymous. These terms refer to a type of life similar to that of the human 
beings, since they are able to develop advanced technologies. 

In this study, we present a chemical evolution model (CEM) for the halo 
and disk components that predicts the temporal behavior of the space distri- 
bution of Z and supernova occurrence in M31, built on precise observational 
constraints (Sec. 2). Based on a set of biogenic, astrophysical, and geophysi- 
cal restrictions, the CEM results lead to the determination of the GHZ (Sec. 

3) . We discuss the implication of the chemical evolution model and the GHZ 
condition on the location, size, and age of the Galactic Habitable Zone (Sec. 

4) . Finally, we present our conclusions of the present study (Sec. 5). 

2. CHEMICAL EVOLUTION MODEL, CEM 

Chemical evolution models study the changes, in space and time, of: a) 
the chemical abundances present in the interstellar medium (ISM) , b) the gas 
mass, and c) the total baryonic mass of galaxies and intergalactic medium. 
These studies have a considerable number of free parameters and observational 
constraints allow us to estimate some of them. If the number of observational 
constraints is high and the observational data are precise, a more solid model 
can be proposed and therefore a better GHZ inferred. In consequence, before 
building the CEM, we collected and determined a reliable chemical data set. 

2.1. Observational constraints 

The present CEM was built to reproduce M31's three main observational 
constraints of the galactic disk: the radial distributions of the total baryonic 
mass, the gas mass, and the oxygen abundance. Since the data come from 
several authors, who have adopted different distances, the data used in this 
work were corrected according to our adopted distance of 783 kpc for M31. 

2.1.1. Radial distribution of the total mass surface density in the disk, Et(^) 

The luminosity profile is produced by the stellar and ionized gaseous com- 
ponents of the galaxies. In evolved spiral galaxies, such as M31, the current 
luminosity is mostly due to the stellar component. Because the luminosity 
profiles of such galaxies follow an exponential behavior with respect to galac- 
tocentric distance (r), we associated that profile to the radial distribution of 
total mass surface density, £y(r) = Y>oexp(— r/ra) where = 5.5 kpc (Renda 
et al. 2005). The value for So is 548.2 M Q pc~ 2 and was obtained by the inte- 
gration of Er(r) over the disk's surface in order to reproduce the total mass 
disk of M31, which is 7.2 x 10 10 M Q (Widrow, Pcrrett & Suyu 2003). 
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2.1.2. Radial distribution of oxygen abundance, [0/Hj(r) 

The data were taken from the compilation of M31's ionized hydrogen (HII) 
region spectra by Blair et al. (1982) and Galarza et al. (1999). From the total 
region sample, we rejected those with uncertain measurements of [Oil], [OIII], 
and [Nil] linesi Taking into account these prescriptions, the final number of 
HII regions of our sample is equal to 83. Once we selected our sample, we 
employed the R23 bi-evaluated method from Pagel et al. (1979), which links 
the intensity of strong [Oil] and [OIII] emission lines with [O /H] , to compute 
the abundances. 

The disadvantage of this method is that it is double-valued with respect 
to metallicity. In fact, at low oxygen abundances ([O/H] < —0.66) the R23 
index decreases with the abundance, while for high oxygen abundances ([O/H] 
> —0.41) the metals cooling efficiency makes R23 drop with rising abundance. 
In order to break the R23 method's degeneracy, we used [Nil] /[Oil] line ratios, 
which are not sensitive to the ionization parameter and are a strong function 
of O/H above log([NII]/[OII]) > 1.2 (Kewley & Dopita, 2002). 

To obtain an estimation of the uncertainties in the abundances, we have 
propagated the error in the line fluxes and added quadratically the accuracy 
of each method, which is between ±0.1-0.2 dex. It is worth mentioning that 
the oxygen abundances derived in this study correspond to the composition of 
the ionized gas phase of the interstellar medium, without taking into account 
oxygen depletion in dust grains. 

In this study, we express the chemical abundances as a function of the 
galactocentric distance (r), related to the Sun as: [0/H](r) = log(0/H)(r) — 
log(0/H)©, where (O/H) represents the ratio of the abundance by number of 
oxygen and of hydrogen, and (O/H)© is this ratio in the Sun, which has a 
value of —3.34 dex (Grevesse, Asplund & Sauval 2007). Moreover, the radial 
abundance gradient is expressed by the slope and y-intercept value of a linear 
relationship that fits the trend of empirical abundances as a function of r. 

The Galactocentric distances of the objects have been derived taking into 
account the Galactocentric distances derived by Blair et al. (1982) and 
Galarza et al. (1999), recomputed for an inclination angle, i = 77°, and 
adopting the distance to M31 from Holland (1998): 783 ± 30 kpc. 

In Figure 1, we compile the radial [O/H] gradients obtained from the dif- 
ferent empirical and theoretical calibrations. It is clear that all the theoretical 
methods give very similar results in contrast to the empirical methods, which 
have much lower y-intercept values even though they keep similar slopes. 

Additionally, we tested the effect of the presence of intrinsic scatter in the 
data in a similar manner to Rosolowsky & Simon's (2008) method for the 
gradient of M33 (the third spiral galaxy of the Local Group) . Following these 
authors' procedure, we applied the method of Akritas & Bershady (1996) to 
compute the gradients in the presence of an intrinsic scatter. We constructed 
histograms for each calibration taking into account samples of ten HII regions 

5 The notation [XI] implies the atom X in neutral form, [XII] indicates X+, [XIII] is X 2+ , 

etc 
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Fig. 1. Gradients obtained using different calibration methods. Empirical methods: 
Pilyugin (2001) and Pilyugin & Thuan (2005). Theoretical methods: McGaugh 
(1991), Kewley & Dopita (2002), Kobulnicky & Kewley (2004), and Zaritsky et al. 
(1994). The most probable gradient among the theoretical methods, applying the 
method of Akritas & Bershady (1996) (see section 2.1.3), is indicated by the broad 
line. That gradient is the one to be reproduced by the chemical evolution model. 
Note that the theoretical methods give very similar results. 

drawn randomly and we obtained the distributions shown in Figure 2. The 
final adopted gradient is the weighted avergare of the four theoretical cali- 
brations we have considered. Hence, the adopted slope is —0.015 ± 0.003 dex 
kpc -1 and the y-intcrccpt value of [O/H] at the center of M31 is 0.44±0.04 dex. 
Finally, we adopted two gradients with which to work: at first, the gradient 
obtained as the most probable value using the method of Akritas & Bershady 
(1996) ([0/H](r) = -0.015 ± 0.003 dex kpc^x r(kpc) + 0.44 ± 0.04 dex), 
which we considered to be representative of the gradients obtained by using 
photoionization model calibrations. On the other hand, we tried to explore 
the chemical evolution of M31 by adopting the gradient given by Pilyugin's 
(2001) empirical calibration, [0/H](r) = -0.014 ± 0.004 dex kpc _1 x r(kpc) 
+ 0.02 ± 0.05 dex), a gradient with similar slope but with a y-intercept value 
0.42 dex lower compared to the most probable gradient obtained from theo- 
retical calibrations. We also tested the gradient given by Pilyugin & Thuan's 
(2005) empirical calibration, which gives a slope 0.006 dex kpc -1 flatter than 
all the other calibrations with an intermediate y-intcrccpt value (see section 
4.1 for a more detailed discussion). 
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Fig. 2. Distribution of the slope (left) and y-intercept (right) of the O/H gradi- 
ents obtained from the four theoretical calibrations considerede in this work: Mc- 
Gaugh (1991) (M91), Kobulnicky & Kewley (2004) (KK04), Kewley & Dopita (2002) 
(KD02) and Zaritsky et al. (1994) (Z94). The distributions are for samples often HII 
regions drawn randomly from the total sample (see text). See figure 5 of Rosolowsky 
& Simon (2008) for a similar diagram in M33. 
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2.1.3. Radial distribution of the gas mass surface density, S gas (r) 

S gas (r) represents all gaseous stages which contain hydrogen, helium, and 
the rest of the chemical elements, i.e. E ffas (r) = + + 'Ez- 

E. M. Bcrkhuijsen (2008, private communication), based on Nieten et al. 
(2006), kindly provided us with the updated atomic and molecular surface 
density of hydrogen for the northern and southern halves of the disk of M31. 
We averaged out these data to obtain Taking into account the He and 
O enrichment by Carigi & Peimbert (2008) (He = 0.25 + 3.3 x O, where He 
and O are abundances by mass), the [O/H] gradient shown in the previous 
subsection, and scaling O and Z to solar values, we computed TiRe and T,z- 
As additional data, we used the compilation of Renda et al. (2005), which 
was corrected following the previously described procedure. 

In Figure 3a we show both sets of corrected data which we used as S gas (j) 
constraints on our chemical evolution model. 

2.2. Model's assumptions 

In the present article, we built a dual-infall model in an inside-out forma- 
tion framework (i.e., the galaxy is formed more efficiently in the inner regions 
than at the periphery), similar to the model used by Renda et al. (2005), 
Hughes et al. (2008), and Carigi & Peimbert (2011) based on the following 
assumptions: 

1. The halo and the disk are artificially projected onto a single disk com- 
ponent of negligible width with azimuthal symmetry; therefore, all func- 
tions depend only on the galactocentric distance (r) and time (t). 

2. M31 was formed by a dual-infall, dT<T(r,t)/dt, of primordial material 
(H = 0.75 and He = 0.25, which are the abundances by mass of hydro- 
gen and helium, respectively) given by the following expression: 

dE T (r,t)/dt = a h (r)e-^ rh + a d (r)e^*- 1G y)/ rd . 

The first term represents the halo formation during the first gigayear 
(Gy), at{r) was obtained considering the halo's present-day total surface 
density profile as 6M Q pc 2 /(l + (?'/8kpc) 2 ) and Th = 0.1 Gy according 
to Renda et al. (2005). 

The second term represents the disk formation from 1 Gy until tg — 13 
Gy (current time), where ad(r) was obtained from £t(0 (see section 
2.1.1), Td = 0.45 (r/kpc) Gy, which represents the inside-out scenario and 
such a value was adopted to reproduce the T, gas (r) along with [0/H](r) 
(see sections 2.1.2 and 2.1.3). 

3. The star formation rate (SFR) is the amount of gas which is transformed 
into stars and was parameterized as the Kennicutt-Schmidt law (Ken- 
nicutt, 1998), SFR(r, t) = vY7^ as (r, t) , where v is the efficiency of star 
formation and n is a number between 1 and 2. During the halo phase, we 
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chose n = 1.0 and we obtained v = 0.50 Gy 1 in order to reproduce the 
maximum average of mctallicity shown by the halo stellar population, 
log(Z/Z ) = -0.5 dex (Koch et al. 2008) (Z Q = 0.012 from Grevesse, 
Asplund & Sauval, 2007). For the disk phase, we chose n = 1.45, a typ- 
ical value for the disk of spiral galaxies (Fuchs, Jahrci/3 & Flynn 2009), 
and we obtained v = 0.23 Gy -1 (M /pc 2 )~ a45 in order to reproduce 
^gas(f) along with [0/H](r) restrictions at the present time. 

4. The initial mass function (IMF) is the mass distribution of the stars 
formed. We considered the IMF of Kroupa et al. (1993), between a 
stellar mass interval of 0.1 to 80 Mq, since the chemical evolution models 
that assume this IMF successfully reproduce the chemical properties of 
the MW's halo-disk (Carigi et al. 2005) . 

5. We adopted the instantaneous recycle approximation, which assumes 
that stars whose initial mass is higher than 1 M die instantly after 
being created and the chemical elements they produce are ejected into 
the interstellar medium. This approximation is good for the elements 
produced mainly by massive stars, such as oxygen, since those stars have 
short lifetimes (10~ 3 to 10~ 2 Gy). The stellar mass fraction returned 
to the ISM by the stellar population and chemical yields are taken from 
Franco & Carigi (2008). 

6. The loss of gas and stars from the galaxy to the intergalactic medium 
was not considered. All chemical evolution models of spiral galaxies 
reproduce the observational constraints without assuming material loss 
from the galactic disk (e.g. Carigi & Peimbert, 2011). Moreover, no 
substantial amount of gas surrounding either M31 or the MW has been 
observed. 

7. We do not include radial flows of gas or stars. 

2.3. Results of the Chemical Evolution Model 

Based on the previously mentioned assumptions and the proper values of 
the free parameters chosen in order to reproduce £t(t*, t g ), Y. gas (r,t g ) and 
the most probable [O/H] gradient (see section 2.1) at the present time, we 
obtained the following results: 

2.3.1. Radial distribution of the gas mass surface density, T, gas 

S gas was solved numerically for each radius at different times. With the 
halo and disk prescriptions shown in section 2.2, the predicted E gas at the 
present time was found and depicted in Figure 3a along with the data. The 
maximum of the theoretical S ffas is shifted towards the inner radius compared 
to the observed values, but the agreement with the central and outer radius 
is quite good. 
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Fig. 3. Results of the chemical evolution model (CEM) at the present-time com- 
pared to the observational data, a) Radial distribution of the gas mass surface 
density (see section 2.1.3). b) [O/H] gradient. The white line is the predicted gra- 
dient; the black line represents the most probable gradient and the shaded area 
compiles the gradients obtained by the theoretical methods (see section 2.1.2). 

We could improve the T, gas agreement for r < 8 kpc, but the slope of 
the chemical gradient (the most important constraint in this study) became 
steeper than the most probable one and it was not precise enough to reproduce 
the chemistry. 

To improve the agreement of E 9a s for r < 8 kpc it may be necessary to 
assume gas flows through the disk towards the galactic center, as an effect of 
the Galactic bar ( Portinari & Chiosi 2000, Spitoni & Matteucci 2011), which 
is beyond the scope of this paper. 

2.3.2. Chemical abundances 

The CEM was built to mainly reproduce the O/H gradient determined 
in this study. That gradient is precise and allowed us to obtain a reliable 
chemical history for M31, where the GHZ will be supported. 

In Figure 3b we present the [O/H] gradient predicted by the model at the 
present time (13 Gy) compared with the observational constraints for M31. 
The gradients obtained by the theoretical methods were enclosed in the shaded 
area of the figure, and the adopted gradient, [O/H] = —0.015 dex kpc^x r 
(kpc) + 0.44 dex, was also included in the figure as a broad white line. The 
gradient obtained from the CEM, [O/H] = -0.015 dex kpc _1 x r(kpc) + 0.45 
dex, represented by a solid black line in the figure, is in perfect agreement 
with the observed value, ensuring the predicted metallicity history. 

As previously described, oxygen is the most abundant elemental compo- 
nent of Z; a chemical evolution model built to reproduce the [O/H] gradient 
therefore allows an adequate approximation of the evolution of the heavy el- 
ements, represented by Z. 
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Fig. 4. Evolution of radial distributions at 1, 3, 5, 7, 9, 11, and 13 Gy. a) Metallicity 
relative to the solar value, log(Z/ZQ) (Zq = 0.012, Grevesse et al. 2007). b) Star 
formation rate, SFR. The horizontal line represents the average value of the SFR 
that the solar neighborhood has undergone during the Sun's lifetime. (Carigi & 
Peimbert 2008). 



Since the GHZ depends mainly on metallicity, we are interested in the 
evolution of the metal radial distribution; hence, in Figure 4a we show the 
behavior of log(Z/Z Q ) with respect to galactocentric distances at different 
times. It is not surprising that the current gradient of log(Z/Z Q ), \og(Z / Zq) 
= —0.015 dex kpc _1 x r (kpc) + 0.41 dex, is similar to the [O/H] gradient 
obtained by the CEM, [O/H] = -0.015 dex kpc _1 X r (kpc) + 0.45 dex, 
confirming that oxygen is the most abundant heavy element and behaves as 
Z. 

From Figure 4a, we notice that the gradient flattens from 3 to 13 Gy, due 
to the inside-out scenario: at the beginning of the evolution the infall was 
relevant in the central regions compared to the outer parts, producing higher 
S gas and SFR, rapidly increasing the oxygen abundance in the inner regions; 
then the infall, and therefore the SFR, dropped, causing a lower increase in the 
oxygen abundance. In contrast, at the final stages of the evolution, the infall 
became important in the outer regions compared to the inner regions, thus 
Ugas increased and SFR became more efficient, producing a higher increase 
in the O abundance in outer regions and, consequently, the gradient became 
flatter. 

The positive slope of the Z gradient at 1 Gy is quite remarkable. This 
opposite slope behavior is due to the enormous amount of primordial material, 
from the intergalactic medium, that fell into the inner parts at the beginning 
of the disk formation, causing gas dilution and a consequentl decrease in the 
O/H values. 
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2.3.3. Supernova rate 

As we have mentioned, the occurrence of supernovae (SNe) is important in 
the calculation of the GHZ because a high SN rate may destroy the atmosphere 
and life of a planet. From the star formation rate inferred by our chemical 
evolution model, wc obtained the SN rate. It was computed according to 
RSN(t) = Nms x SFR(i), where N MS = 0.05MQ 1 is the number of type II 
supernova progenitors per Mq. We assumed that those progenitors are stars 
more massive than 8 Mq and therefore Nms was calculated by integrating 
the initial mass function between 8 and 80 Mq. 

In Figure 4b, we show the star formation rate as a function of galactocen- 
tric distance for the same times shown in Figure 4a. In order to compare with 
the SFR that occurred in the solar neighborhood during the Sun's lifetime, 
we added the average value of the SFR in the solar neighborhood over the last 
4.5 Gy. That average was obtained from the mean SFR at the solar radius 
during the last 4.5 Gy, <SFR (8 kpc)>=<SFR sy >= 3.8 M Q Gy" 1 pc" 2 
(see Fig. 2 by Carigi & Peimbert 2008). Consequently, the 4.5 Gy average 
SN rate of the solar neighborhood is < RSN S v >= 0.2G?/™ 1 pc~ 2 . 

3. GALACTIC HABITABLE ZONE 

The Galactic Habitable Zone (GHZ), is defined as the region with sufficient 
abundance of chemical elements to form planetary systems in which Earth-like 
planets could be found and might be capable of sustaining life. Therefore, a 
minimum metallicity is needed for planetary formation, which would include 
the formation of a planet with Earth-like characteristics (Gonzalez et al. 2001, 
Linewcaver 2001), and a CEM provided the evolution of the heavy element 
distribution (see previous section). 

3.1. Characteristics and conditions for the GHZ 

In this context, it is important to bear in mind the properties and char- 
acteristics that an Earth-like planet and its environment need in order to be 
considered habitable or leading to habitability. In astrobiology, the chemical 
elements are classified according to the roles they play in the formation of 
an Earth-like planet and in the origin of life: the major biogenic elements 
are those that form amino acids and proteins (e.g. H, N, C, O, P, S), the 
geophysical elements are those that form the crust, mantle, and core (e.g. Si, 
Mg, Fe) (Hazcn ct al. 2002). 

3.1.1. Astrophysical conditions 

The conditions of the environment around an Earth-like planet are the 
following: 
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1. Planet formation: The metallicity in the medium where the planet may 
form (protoplanetary disk) must be such that it allows matter conden- 
sation and hence protoplanet formation (Lineweaver 2001). Therefore, 
Earth-like and Jupiter-like planets may be created. Whether the cre- 
ation of Jupiter-like planets is beneficial or not is a debatable subject 
(Horner & Jones 2008, 2009). They may act as a shield against mete- 
oritic impact (Fogg & Nelson 2007, Ward & Brownlee 2000), or might 
migrate towards the central star of the planetary system, thereby af- 
fecting the internal planets, or even destroying them (Lineweaver 2001, 
Lineweaver et al. 2004). On the other hand, the simulation by Ray- 
mond ct al. (2006) predicts that Earth-like planets may form from 
surviving material outside the giant planet's orbit, often in the hab- 
itable zone and with low orbital eccentricities. Most (~ 90%) of the 
planet-harboring stars, detected by several research teams, show a wide 
range of stellar metallicity, —0.35 < \og(Z/Zo) < +0.35, and the dis- 
tribution peaks at \og{Z / Zq) ~ +0.20 (See Figure 5). Data compila- 
tion taken from The Extrasolar Planets Encyclopaedia on July 2012, 
http : / /exoplanet . eu / catalog. php[ ) . Moreover, based on the same data 
source, the observed relation between stellar Z abundance and the planet 
mass presents a wide dispersion. Since that sample includes all types of 
planets — both Earth-like and Jupiter-like — the stellar Z— planet mass 
relation does not show a clear trend. Recently, Jenkins ct al. (2012) 
have grown the population of low-mass planets around metal-rich stars, 
increasing the dispersion in the stellar Z— planet mass relation. There- 
fore we think that the log(Z/Z©) range represents the necessary Z to 
form planets that survive to planetary migrations. It should be men- 
tioned that the stellar Z is associated with the planet Z since the planets 
were created from the same gaseous nebula of the star. 



2. Survival from SNc: When a supernova explodes, it emits radiation capa- 
ble of wiping out the atmosphere of a nearby planet and sterilizing the 
life that may lie on its surface (Lineweaver et al. 2004, Gowanlock et al. 
2011). Since the Earth is the only known planet with life, we assumed 
it represents the survival pattern to SN explosions. We considered a 
planet to survive if the rate of SN (RSN) exploding, within a distance 
of 0.1 kpc, is the same or less than the time-average, during the last 4.5 
Gy of SN rate in the solar neighborhood (< RSNsv >), where the Sun 
and its planets formed. That number is < RSNsv >= 0.2 Gy~ pc -2 
(see section 2.3.3). We have analyzed two possibilities for the planet 
survival: i) the life on the planet will be annihilated forever when the 
instantaneous RSN(t) is higher than < RSNsv >, and ii) the exoplanet 
will be sterilized when the average RSN, during the entire existence of 
the planet, is higher than < RSNsv >■ 
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Fig. 5. The fraction of detected exoplanets as a function of the metallicity of 
the stars with exoplanets. Data compilation taken from The Extrasolar Planets 



Encyclopaedia on July 2012, http://exoplanet.eu/catalog.php 



3.1.2. Geophysical characteristics 

An Earth-like planet is one that has enough biogenic and geophysical el- 
ements to sustain and allow the development of life as it is known on Earth 
(Sleep, Bird & Pope 2012, Bada 2004, Gomez-Caballero & Pantoja-Alor 2003, 
Hazen et al. 2002, Orgel 1998). In general terms, for a planet to be identified 
as an Earth-like planet, it has to satisfy the following conditions: 

1. Tectonic plates: a crust formed mainly of Si constitutes the tectonic 
plates. These are necessary since life must have a site to live on with 
enough resources for survival. The recycling of the tectonic plates keeps 
the density and the temperature of the planetary atmosphere and the 
amount of liquid water and carbon needed for life. (Hazen et al. 2002, 
Lineweaver 2001, Scgura & Kaltenegger 2006 ). 



2. Water: on the surface of the tectonic plates are found the oceans, formed 
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by H2O. Since life on Earth is the pattern for life, water is a crucial 
resource for the emergence and survival of life (McClcndon 1999). 

3. Atmosphere: the atmosphere of a planet should be dense enough to 
protect the planet from UV radiation and meteoric impacts, and thin 
enough to allow the evolution of life on the planet's surface. The atmo- 
spheric composition of the primitive Earth allowed the origin of life and 
had abundance of CO, CO2, H 2 0, N 2 0, and N0 2 . Such compounds 
were crucial to the origin of the present atmosphere on Earth. (Bada 
2004, Chyba & Sagan 1991, Maurette et al. 1995, Navarro-Gonzalez et 
al. 2001, Sekine et al. 2003). 

Since all the geophysical elements are heavier than He, the abundance 
of these kinds of elements necessary to create an Earth-like planet is both 
contained and well represented by the distribution of Z given by the chemical 
evolution model. 

3.1.3. Biogenic characteristics 

As previously mentioned, the GHZ is based on the pattern of life on Earth; 
therefore, in order to determine the time of the origin and development of life 
on Earth several studies were taken into account. 

There are several theories of the means and places where life originated on 
Earth. Some of those theories are based on hydrothermal vents on the seabed, 
where an interaction between the terrestrial mantle and the ocean exists (Bada 
2004, Chang 1982, Gomez-Caballero & Pantoja-Alor 2003, Hazen et al. 2002, 
McClcndon 1999). Other theories have assumed the migration of life from 
other regions of space to Earth (Maurette et al. 1995, Orgel 1998). There 
are others explaining the catalytic effect of lightning or metallic meteorites 
(Chyba & Sagan 1991, Navarro-Gonzalez et al. 2001, Sekine et al. 2003) in 
the primitive atmosphere and/or in the seas. Regardless of how life emerged, 
all theories have agreed on the same biogenic requirements, i.e. C, N, O, P, 
S. Those biogenic elements are heavier than He and the abundance of these 
elements is consequently contained and well represented by the distribution 
of Z given by the CEM. 

The earliest evidence we have for life on Earth is about 3.5 Gy ago (the 
most ancient fossils being cyanobacteria) or 3.8 Gy (contested carbon isotope 
evidence). There is general agreement that life got started 3.7 Gy ago (Ricardo 
& Szostak 2009). The the Earth is 4.6 ± 0.1 Gy old (Bonanno et al. 2002); 
therefore, life took around 0.9 Gy to emerge. In consequence, in this study 
we consider 1.0 Gy to be the minimum age of a planet in order for basic life 
to appear on its surface. 

On the other hand, if evolved life is associated with humans, another 
parameter of life should be taken into account since humans appeared on 
Earth 2 million years (2 x 10~ 3 Gy) ago (Bada 2004). However, the fact that 
the times used in our CE study arc of the order of gigayears, this 2 x 10~ 3 
Gy is negligible. Thus, in this study we consider 4.5 Gy to be the minimum 
age of a planet for complex life to evolve. 



CHEMICAL EVOLUTION AND GHZ OF M31 



15 



3.2. Restrictions for the GHZ 

Based on the previously collected information, in order to obtain the evo- 
lution of the Galactic Habitable Zone in M31, we chose the following astro- 
nomical and biogenic restrictions: 

1. Regions in the galaxy with metallicity between —0.35 < \og{Z / Zq) < 
+0.35 are able to form most Earth-like planets. 

2. Earth-like planets require 1.0 Gy to create basic life (BL). 

3. Earth-like planets need 4.5 Gy to evolve complex life (CL). 

4. Life on formed planets is annihilated by the supernova explosions under 
the following conditions: 

(a) The SN rate at any time has been higher than the average SN rate 
in the solar neighborhood (< RSNsv >) during the last 4.5 Gy of 
the Milky Way's life, and 

(b) The average RSN, during the whole existence of the planet, is 
higher than < RSN S v > 

Based on: i) Planet formation, dependent on Z, is the crucial factor of the 
GHZ (Suthar & McKay 2012) and ii) the entire galactic disk has undergone 
considerable star formation, where and when the Z values are suitable for 
planetary formation (see Figures 4a, 4b, and 5), then the probability, P(r,t), 
to form Earth-like planets which survive SN explosions and where basic and 
complex life may surge is calculated as 

P(r,t) = P z (r,t) x P BL (t) x P CL (t) x P 8N {r,t), where 

1. Pz(r,t) is the Z-depcndent probability of forming terrestrial planets. 
We adopt a Z dependence for Pz(r,t) identical to the Z distribution 
shown by extrasolar planets (see Figure 5). 

2- Pb l {t) and PcLif) are the probabilities of the emergence of basic life 
and the evolution of complex life. We adopt functions with no dispersion 
(equal to 0.0 or 1.0) because the uncertainties in the evidence for life 
are similar to the temporal resolution of the chemical evolution model. 

3- PsN{r,t) is the probability of survival of supernova explosions. We 
adopt PsN{r,t) = 1.0 or Psiv(r,t) = 0.0 if SN rate is lower or higher 
than the average RSN, respectively. This is a good approximation 
because of the quick decay of Psn (Lineweaver et al. 2004) and the two 
extreme conditions supported on the pattern of life and considered in 
this study. 
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3.3. Results of the Galactic Habitable Zone 

Our chemical evolution model — explained in section 2 — provided the evo- 
lution of metals and the SN rate at different galactocentric distances, which 
are the data required to compute the Galactic Habitable Zone. We then ap- 
plied the restrictions listed in the previous sections on the galactic disk of 
M31 and obtained the GHZ. Since most previous studies of the GHZ have fo- 
cused on the Milky Way (Lineweaver et al. 2004, Prantzos 2008, Gowanlock 
et al. 2011), we also obtained the GHZ of the Andromeda galaxy using their 
common restrictions in order to compare the GHZ in both neighboring spiral 
galaxies using identical galactic habitable conditions. 

3.3.1. Present work conditions 

Based on the results of the chemical evolution model, we found the times 
and regions in the galactic disk where the metallicity condition is satisfied; 
we then applied the condition of supernova survival, and imposed the restric- 
tion of basic life (1 Gy) and evolved life (4.5 Gy). In Figure 6 we show the 
evolution of the GHZ assuming these constraints and we mark with differ- 
ent symbols the zone where there are enough metals to form planets (green 
area), where planets were sterilized (red), where planets may have basic life 
(striped pattern), and where planets might harbor complex life (cross-hatch 
pattern). We indicate with differently shaded levels the probability values of 
forming Earth-like planets. The yellow and blue areas represent the zones 
with probabilities lower than 4%, because the metallicity is too high and too 
low, respectively. Terrestrial planets could hardly be formed in these zones. 

The only difference between the models shown in Figure 6a and Figure 
6b is the assumption about planet survival with respect to the occurrence of 
SNe. In Figure 6a we assumed that life survives if RSN(t) is smaller than 
< RSNsv > at each instant of the planet's existence. This implies that life 
on the planet would not recover after a nearby SN (~ 6 supernovae per million 
year). In Figure 6b, we considered life to survive if the average SN rate during 
the entire existence of the planet is lower than < RSNsv >• This implies 
that life might recover or reappear after a nearby SN explosion. 

The GHZs shown in Figure 6 are smaller at high planet ages (equivalent 
to low evolutionary time) and excludes short galactocentric radii, due to the 
inside-out formation scenario, upon which the chemical evolution model was 
built. That model predicts a high star formation rate at inner radii during the 
first moments of evolution and, consequently, a high SN rate that would have 
sterilized the planets formed there. If life has not recovered after RSN(t) =< 
RSNsv >, then there can be no planets with complex life at r < 13 kpc (see 
Figure 6a), but if life does recover, there could be some planets with complex 
life at any r (see Figure 6b). 

Based on both figures, the GHZ of M31 is located in the galactic disk and 
consists of a ring that widens from the central to the outer parts, from the first 
Gy to 5 Gy with a low probability of forming terrestrial planets. During the 
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Fig. 6. Evolution of the Galactic Habitable Zone of M31 assuming our considera- 
tions (see section 3.2). Green area: zone with sufficient metallicity to form planets 
according to the Z distribution shown by detected stars with extrasolar planets. The 
width at half maximum of this distribution is located at —0.35 < log(Z/Z0) < +0.35 
and the peak is at \og(Z/Z@) ~ 0.2. Shaded areas: probability values of forming 
Earth-like planets. Yellow and blue areas: zones with probabilities lower than 4% 
due to the metalicity being too high and too low, respectively. In these zones terres- 
trial planets have hardly formed. Striped pattern: zone where planets may harbor 
basic life. Cross-hatched pattern: zone where planets may harbor evolved life. Red 
area: zone where life may be extinguished by SN explosions, comparing two dif- 
ferent regimes of SN rate to the time-average SN rate of the solar neighborhood 
during the Sun's existence, < RSNsv >■ Left: if planet sterilization occurs when 
the instantaneous SN rate is higher than < RSNsv >■ Right: if planet sterilization 
occurs when the average SN rate during the planet's entire existence is higher than 
< RSNsv >■ 



next ~ 2 Gy, the ring expands into the inner parts, reaching practically the 
whole galactic disk, between 3 and 30 kpc, increasing the probability towards 
smaller radii. After 7 Gy, the high probability shifts to central radii. 

The ring contains planets not older than 11 Gy, those planets with ages 
higher than 1 Gy harbor basic life, but evolved life might be found on planets 
older than 4.5 Gy. If life is annihilated with an instantaneous SN rate higher 
than what the Earth has experienced, most planets formed between 3 and 13 
kpc and older than 2 Gy would have been sterilized, and most the probable 
GHZ with complex life would be located between 13 and 24 kpc. If life is 
extinguished with an average SN rate over the entire existence of the planet 
higher than that undergone by the Earth, planets formed between 3 and 10 
kpc and older than 6.5 Gy would be sterilized, and the most probable GHZ 
with complex life would be located between 3 and 13 kpc. In either case, the 
halo component (for evolution times lower than 1 Gy at any r) is discarded 
from the GHZ because of its low Z. 
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3.3.2. Milky Way conditions 

Four of the five previous studies of the GHZ (Gonzalez et al. 2001, 
Line-weaver et al. 2004, Prantzos 2008, Gowanlock et al. 2011, Suthar & 
McKay 2012) have focused on the Milky Way. The Milky Way is a smaller 
spiral galaxy than the Andromeda galaxy, with a lower halo metallicity, a 
steeper slope, and a lower y-intercept value of the O/H gradient; hence, the 
MW has had a chemical history different from that of M31. The compari- 
son between the GHZs of the MW and M31, from the same astrophysical, 
geophysical and biogenic restrictions, have allowed us to explore the effect on 
the GHZ of the intrinsic differences among spiral galaxies. These differences 
are mainly the infall and the star formation history, which determine the gas 
enrichment in each galaxy. 

Lineweaver et al. (2004), Prantzos (2008) and Gowanlock et al. (2011) 
have computed the GHZ of the MW assuming specific Z ranges to form Earth- 
like planets and for SN survival. Since Gowanlock et al. (2011) adopted a 
much wider Z range than that of Lineweaver et al. (2004) and with probabil- 
ities lower than that considered by Prantzos (2008), we computed the GHZ 
of M31 based on the assumptions of Lineweaver et al. (2004) and Prantzos 
(2008), separately. 

First, we computed the GHZ of M31 using Lineweaver et al.'s assumptions: 

i) Earth-like planets, that survive the migration of Jupiter-like planets, 
formed in regions where —0.35 < log(Z/Z Q ) < +0.35, a range that represents 
the width at half maximum of the distribution considered by Lineweaver et 
al. and peaks at \og(Z / Zq) ~ 0.1. This Z distribution is poorly populated at 
high metallicities, according to extrasolar planets data up till 2004. 

ii) Those planets formed are sterilized by the SN explosions if the time- 
average of the neighboring SN rate during the first 4.5 Gy of the planet's 
age, < RSN(r,t) >, is higher than twice that of the mean SN rate in the 
solar neighborhood during the Earth's existence, < RSNsv >■ Specifically, 
PsN{r 7 t) = 1.0 or PsN(r,t) = 0.0 if the average SN rate, during the first 4.0 
Gy of the planet's existence, is lower or higher than 2x < RSNsv >■ 

In Figure 7a we show the M31 GHZ, taking into account the previous 
conditions, and also we impose our basic and complex life conditions. This 
GHZ is similar, but not equal, to that computed with our conditions (see 
Figures 6) because its Z range is identical to ours, but its probability of 
forming Earth- like planets and its survival condition are slightly different. The 
central area of the galactic disk, discarded by planet sterilization, is wider than 
that in Figure 6b, because its survival condition (average in the first 4.0 Gy) 
implies that life is less resistant to SN explosions than our survival condition 
(averaged over the planet's entire existence). That fact, together with the low 
probability of forming planets at high Z, places the most probable GHZ with 
complex life between 10 and 19 kpc. 

In Figure 7b we show the M31 GHZ taking into account the Z and survival 
assumptions by Prantzos (2008) plus our basic and complex life conditions. 
The only difference with the Lineweaver et al. conditions is the Z distribution 
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Fig. 7. Evolution of the Galactic Habitable Zone of M31 assuming the GHZ consid- 
erations for the MW (see section 3.3.2). Green area: zone with sufficient metallicity 
to form Earth-like planets that survive to migration of Jupiter-like planets. Shaded 
areas: probability values of forming those planets. Yellow and blue areas, and striped 
and cross-hatched patterns as in Figure 6. Red area: zone where life may be extin- 
guished if the average SN occurrence, during the first 4.0 Gy of the planet existence, 
is higher than 2x < RSNsv >■ Left: Lineweaver et al. (2004) ^-distribution condi- 
tion: -0.35 < log(Z/Z Q ) < +0.35 with a maxi mun at log(Z/Z ) = +0.10. Right: 
Prantzos (2008) ^-distribution: -1.00 < log(Z/Z Q ) < +0.50 with equal probability 
for each Z. 

and range to form Earth-like planets that survive to the migration of Jupiter- 
like planets. Prantzos assumed a wider range in Z, —1.00 < log(Z/Z Q ) < 
+0.50, with an almost constant probability for any Z. Consequently, the 
habitable galactic zone includes the entire galactic disk since the first Gy with 
equal probability for each radius and time. Earth-like planets contained in this 
zone are as old as 12.5 Gy. In this case, the GHZ extends to the innermost 
regions of the disk, where Z is the highest, and to the youngest halo star, 
where log(Z/Z Q ) 1. 

4. DISCUSSION 
4.1. Chemical Evolution model 

A chemical evolution model can be computed with three minimum obser- 
vational constraints at the present time: total baryonic mass, gas mass, and 
the abundance of one chemical clement. In this study we have compiled from 
the updated literature and well described T>t(t) and E gas (r) distributions, 
and we have determined a precise O/H gradient from HII regions. Based 
on these quite restrictive constraints, we obtained a solid chemical evolution 
model of the halo and disk of M31. 

The chemical evolution model was built to match the most probable O/H 
gradient, [0/H](r) = —0.015 dex kpc" 1 x r(kpc) + 0.44 dex, determined con- 
sidering the intrinsic scatter and several theoretical calibrations (see section 
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2.1.2). This gradient is caused by the inside-out scenario: the infall is faster 
and more abundant in inner regions, producing a more efficient SFR, which 
increases the O/H abundance of the inner parts faster than that of the outer 
parts. 

In this study, we also explored the gradient determined by empirical cali- 
brations and proposed by Pilyugin (2001) and Pilyugin & Thuan (2005). The 
gradient obtained through the Pilyugin (2001) calibration presents the same 
slope, but its y-intercept value is lower by 0.42 dex than that obtained from 
theoretical calibrations. On the other hand, the gradient obtained by using 
the Pilyugin & Thuan (2005) calibration, [0/H](r) =-0.008 dex kpc -1 x 
r(kpc) — 0.23 dex, is 0.007 dex kpc -1 flatter and the y-intercept value lower 
by 0.67 dex than that obtained from theoretical calibrations. We were not 
able to reproduce the Pilyugin gradient without overestimating the Y> gas (r) 
distribution. 

Pcimbert et al. (2007) suggested that the abundances derived from the R23 
empirical calibrations could be underestimated due to the presence of tem- 
perature variations in the ionized volume, which leads to the overestimation 
of the electron temperatures, and that the use of empirical calibrations based 
on optical recombination line measurements (which are much less affected by 
temperature variations and give abundances similar to those obtained using 
theoretical calibrations) could reconcile the R23 theoretical calibrations with 
empirical ones. This is a well known problem in the astrophysics of pho- 
toionized nebulae (see Peimbert et al. 2007 and references therein for more 
discussion about this topic). 

The slope of the O/H gradients which we computed are in some way 
steeper, but are consistent, within the errors, to those derived from super- 
giants within 12 kpc of the galactic center reported by Trundle et al. (2002), 
which amounts to -0.006± 0.020 dex kpc -1 . Venn et al. (2000) had also 
found that there was no oxygen gradient between 10 kpc and 20 kpc; how- 
ever, their study was based on abundance determinations for only three A-F 
supergiants and their results are therefore inconclusive. Kwitter et al. (2012) 
have recently computed a radial [O /H] gradient from planetary nebulae (PNe) 
data. They determined a gradient of — 0.011±0.004 dex kpc -1 , which is con- 
sistent to within the errors with our value derived from H II regions. 

We need additionally to bear in mind that, although the values of the gra- 
dient we have obtained are somewhat shallower than typical values for spiral 
galaxies, recent results seem to confirm that M31 is a barred spiral galaxy 
(see Beaton et al. 2007, Athanassoula & Beaton 2006 and references therein), 
which produce shallower gradients than normal spiral galaxies (Friedli et al. 
1994). 

The S g as obtained from the chemical evolution model follows a similar 
pattern to that observed at radii larger than 9 kpc. However, at smaller 
radii the model predicted more gas than observed. Since T, gas , at a fixed 
galactocentric distance, decreases due to the SFR and increases mainly due to 
infall, we tried to reproduce at smaller r values the T, gas increasing the SFR, 



CHEMICAL EVOLUTION AND GHZ OF M31 



21 



but the slope of the O/H gradient became steeper than the most probable 
one. In addition, we tried to solve the discrepancy between the model and the 
observed T, gas at low radii, changing the radial distribution of the infall, but 
the agreement with the observed E-r was lost. It is possible to reduce S 9as in 
the inner part by adopting gas inflows to the galactic bar ( Portinari & Chiosi 
2000, Spitoni & Matteucci 2011). That assumption will be considered in a 
future article. 

The predicted distribution of the SFR at 13 Gy agrees with the observed 
one (see data collected by Yin et al. 2009 and Marcon-Uchida et al. 2010) 
for inner regions, but is in disagreement for outer regions. On the other 
hand, the agreement between the predicted and observed T, gas is good for 
the outer regions, but not for inner regions. Since the SFR assumed in this 
study depends on the gas mass, there is an inconsistency between the observed 
Ti gas and the SFR. Therefore, it is not obvious how to change the two free 
parameters of our SFR (see section 2.2) to match both the T, ga8 and the SFR 
distributions with the observed values. In order to improve the agreement 
of the obtained Y, gas and SFR, we needed to consider different accretion and 
star formation laws based on more complex physical properties for galactic 
and star formation. 

Renda et al. (2005), Matsson (2008), Yin et al. (2009), and Marcon- 
Uchida et al. (2010) have computed CEMs for the Andromeda galaxy taking 
into consideration different assumptions and using diverse codes. These mod- 
els were built to reproduce different O/H gradients, which show higher dis- 
persion or error, with O/H differences at a given r value in the 0.4 to 0.7 dex 
range; consequently, these models are not constrained as well as our model. 
In particular the agreements of the Marcon-Uchida et al. model are poorer or 
marginal. 

We also made a comparative study of the chemical evolution models and 
the observational constraints of the Milky Way galaxy with the Andromeda 
galaxy, and we noticed that: 

1. The maximum average of the Fe/H abundances present in the halo stars 
of M31 corresponds to log(Z/Z©)= —0.5, while in the MW it amounts 
to —1.6 (Koch et al. 2008). In order to reproduce the higher Z value 
in the halo stars of M31, our model requires, in the halo phase, a more 
efficient infall (r^ is about 5 times lower) and SFR (two orders of mag- 
nitudes higher) than those for the MW model computed with similar 
assumptions (Carigi et al. 2005, Carigi & Peimbert 2011). 

2. The adopted gradient for M31, [0/H](r) = -0.015 dex kpc" 1 x r (kpc) 
+ 0.44 dex, is 0.03 dex flatter and 0.06 dex higher in the central value 
than the O/H gradient of the Milky Way disk determined by Esteban 
et al. (2005) from optical oxygen recombination lines. In order to re- 
produce the flatter M31 gradient, our model needed, in the disk phase, 
a less marked inside-out formation, i.e. slower infall in the inner regions 
and faster accretion in the outer parts compared to the infall considered 
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in CEMs of the MW. Moreover, to reach the higher O/H values, the 
model required a more active star formation history than that assumed 
for the MW disk (Carigi & Peimbert 2011), resulting in a lower T, gas (r) 
in M31 than in the MW, in agreement with the observed distribution in 
both galaxies. 

4.2. Galactic Habitable Zone 

The chemical evolution model discussed above provides the evolution and 
location of the heavy elements of which planets are formed and the rate of 
neighboring SNe that may sterilize life on the planets. The age and the 
location of Earth-like planets inside a galaxy depend mainly on the Z range 
assumed to form those planets. 

It is difficult to quantify the role of metallicity in planet formation. At 
present, theoretical models fail to explain in detail the physical processes in 
protoplanetary disks and in particular the metallicity effect on the planet mass 
(rocky like Earth or gaseous like Jupiter). Furthermore, the observational data 
of detected extrasolar planets do not show a clear trend between the stellar 
metallicity and the planet mass that orbit around it. Recent studies (e.g. 
Jenkins et al. 2012) have increased the dispersion in the stellar Z— planet mass 
relation to detect more low-mass planets around metal-rich stars. Also, it is 
noticeable that in many detected extrasolar planetary systems, Earth-like and 
Jupiter-like planets coexist (The Extrasolar Planets Encyclopedia). Moreover, 
Raymond et al. (2006) have simulated terrestrial planet growth during and 
after giant planet migration. Observational studies on the circumstellar disk 
in young stellar clusters of low and solar metallicities suggest that the disk 
lifetime shortens with decreasing thereby reducing the time available for 
planet formation (Yasui et al. 2010). Based on circumstellar disk models and 
dust grain properties, Johnson & Li (2012) conclude that the first Earth-like 
planets likely formed from gas with metallicities Z > 0.1 Zq. 

For these reasons we adopted the Z distribution shown by planet-harboring 
stars as the metallicity condition for the GHZ instead of undertaking a prob- 
ability study of planet formation and migration as a function of metallicity, 
as did Lineweaver (2001), Prantzos (2008), and Gowanlock et al. (2011). 
Using that Z distribution, the GHZ currently extends in the galactic disk 
from 3 to 30 kpc with maximum planet ages of 11 and 8.5 Gy, respectively. 
The most probable GHZ, considering only basic life, is between 6 and 17 kpc 
and ages between 1 and 4.5 Gy; assuming complex life, however, the most 
probable GHZ is between 3 and 13 kpc and ages between 4.5 and 6.5 Gy 
(see Figure 6b). When we adopted the Z distribution by Lineweaver, with 
a lower number of exoplanets at high metallicities, the most probable GHZ, 
considering only basic life, is between 13 and 24 kpc and ages between 1 and 
4.5 Gy, but, assuming complex life, the most probable GHZ is between 10 
and 19 kpc and ages between 4.5 and 8 Gy (see Figure 7a). In both cases, 
the formation of Earth-like planets was discarded in the Andromeda halo 
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since that galactic component contributes with very low amounts of metals 
(\og(Z ha i /Z Q ) < -0.5, see Figure 4a). 

When we adopted the Z condition by Prantzos, ( — 1.00 < log(Z/Z Q ) < 
+0.50), the first planets formed are older than the galactic disk and a few 
Earth-like planets might even be formed in the halo because the maximum 
Z shown by the halo stellar population is 0.5 dex higher than the minimum 
Z required (\og(Z / Zq) = —1.00) to form Earth-like planets. It is important 
to remark that only in this case have Earth-like planets surviving Jupiter-like 
ones been formed in the innermost regions, because the maximum mctallicity 
reached by the model is log(Z/Z ) = +0.4 (see Figure 7b). 

The age and the location of Earth-like planets sterilized by SN explosion 
inside a galaxy is determined by the assumed SN condition. That condition is 
the most uncertain restriction, because it is difficult to estimate the resistance 
of life to SN explosions, which depends on unknown biogenic and astrophysical 
factors, such as the minimum SN-planet distance that allows life to survive, 
the type of atmosphere and oceans where life could evolve, and the biological 
properties of the different types of life, among others factors. 

Since our life pattern is based on one Earth, life that has been able to 
recover from neighbor SN events, and the solar neighborhood has undergone 
a nearby constant and low SN rate from its formation (Renda et al. 2005, 
Prantzos 2008, Carigi & Peimbcrt 2008), we have assumed that the average 
SN rate of the solar neighborhood is the maximum one that life can with- 
stand during the whole planet evolution. Using that SN restriction, planets 
older than 6.5 Gy and located at r <10 kpc were sterilized (see Figure 6b). 
When we assumed the SN condition by Lineweaver et al. (2004), a less re- 
strictive constraint that considers the time-average of the SN rate during the 
first 4.0 Gy life of the planet, only planets older than 6 Gy located at r<13 
kpc have been sterilized (see Figure 7a). In all cases, even though the area 
subtracted from the GHZ due to SN explosion is small, the number of planets 
where life has disappeared is significant because the stellar mass decreases 
with the galactocentric radius, following an exponential behavior given by the 
luminosity profile (see section 2.1.1). 

In a future paper, the SFR(r,t) will be considered to compute the GHZ 
as the most likely zone to have stars with habitable planets. In that case, 
the GHZ will be located even further towards the centre of the disk of M31, 
similarly to the GHZ of MW obtained by Prantzos (2008) and Gowanlock et 
al. (2012). Since the stellar mass profile of M31 is flatter than that of the 
MW, it will be interesting to know how wide the M31 GHZ is and to quantify 
the number of planet-harboring stars. 

Using the same astrophysical, geophysical, and biogenic restrictions in 
M31 as those considered in the MW studies (Lineweaver et al. 2004, Prantzos 
2008), we found that the GHZ of M31 is wider compared to that of the MW. 
This result is mainly due to the higher and flatter Z gradient of the disk of 
M31. Consequently, the outer parts of the disk are more metal- rich and able 
to form planets. In addition, we found that the GHZ of M31 is older due to 
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the higher metallicity of the halo, causing the first stellar generations of the 
disk to have more heavy elements to form planets than the first stars of the 
MW (See Figures 6). 

Throughout this section we have emphasized that the chemical content is 
the most restrictive condition to the GHZ. Furthermore, we have mentioned 
that the slope of the Z gradient is determined by the efficiency of inside-out 
scenario, whereas the y-intercept values of the gradient depend on by the 
amount of accreted gas and the star formation efficiency. Therefore, we stress 
that the manner in which a galaxy is assembled and its stars formed determine 
the size and age of the Galactic Habitable Zone. 

5. CONCLUSIONS 

Based on O/H values of HII regions, we obtained O/H gradients with 
different empirical and theoretical calibrations. In the presence of intrinsic 
scatter, we computed the most probable gradient for M31 from theoretical 
calibrations based on the R23 method and we found that [O/H] (7-) = — 0.015 ± 
0.003dex kpc- 1 x r(kpc) + 0.44 ± 0.04 dcx. The slope of the gradient is 0.03 
dex kpc -1 flatter, and the value at the center of the galaxy is 0.06 dex higher, 
than those values of the O/H gradient of the Milky Way disk. 

The chemical evolution model built to reproduce our precise O /H gradient 
of the galactic disk matches the radial distribution of the gas mass, S gQS (r), 
for the outer regions quite well, but it fails for the inner ones. On the other 
hand, the radial distribution of the star formation rate, SFR(r), fits for the 
inner parts, but not for the outer ones. Therefore, in order to improve the 
agreements, the model would require more complex galactic and star forma- 
tion histories. 

The model cannot reproduce the O/H gradient computed by empirical 
methods, which is 0.42 dex lower for the central value than the gradient ob- 
tained by theoretical methods, unless the other observational constraints are 
considerably modified. 

Based on the chemical evolution model we obtained, for the first time, the 
Galactic Habitable Zone (GHZ) of M31, considering two space requirements: 
regions with i) sufficient metallicity for planet formation with a probability 
law that follows the Z distribution shown by exoplanets and ii) an average 
SN rate similar to that permitting the existence of life on the Earth; and two 
time requirements: the existence of basic life (like cyanobactcria) , and the 
development of evolved life (like humans). 

The GHZ of M31 is located in the galactic disk and consists of a ring 
that widens from the inner to the outer parts, during the last ~11 Gy of 
evolution, reaching at the present time practically the whole galactic disk. If 
life is extinguished with an average SN rate higher than that for the Earth 
during the planet's entire existence, planets formed between 3 and 10 kpc and 
older than 6.5 Gy would have been sterilized. If only basic life is considered, 
the most probable GHZ corresponds to 28% of the galactic disk area and this 
GHZ is located between 6 and 17 kpc on planets with ages between 4.5 and 1 
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Gy. If complex life is considered, the most probable GHZ corresponds to 18% 
of the galactic disk area and is located between 3 and 13 kpc on planets with 
ages between 6.5 and 4.5 Gy. 

In GHZ studies, the Z restriction is crucial for finding the location of 
Earth-like planets, but the SN survival condition allows us to compute the 
location of those planets that survive SN events. 

Since M31 is more metal-rich than the MW (specifically the metallicity 
of the halo is higher and the chemical gradient of the galactic disk is flatter 
and higher compared to the Galaxy), the GHZ of M31 is wider and older than 
that of the MW. Therefore, the number of planets harboring life in Andromeda 
could be higher than in our galaxy. 

Part of this work was submitted by Sofia Meneses Goytia to the Master's 
Programme in Sciences Chemistry at the Universidad National Autonoma de 
Mexico. 
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